arXivil509.04124v2 [nucl-th] 17 Sep 2015 


Two-particle correlations in psendorapidity in a hydrodynamic model 

Piotr Bozek/’[^ Wojciech Broniowski,^’and Adam Olszewski^ ’[^ 

^AGH University of Science and Technology, Faculty of Physics and 
Applied Computer Science, al. Mickiewicza 30, 30-059 Krakow, Poland 
^ The H. Niewodniczahski Institute of Nuclear Physics, 

Polish Academy of Seiences, 31-342 Krakow, Poland 
^Institute of Physics, Jan Kochanowski University, 25-406 Kielce, Poland 
(Dated: 14 September 2015) 

Two-particle pseudorapidity correlations of hadrons produced in Pb-|-Pb collisions at 
^snn = 2.76 TeV at the CERN Large Hadron Collider are analyzed in the framework of a model 
based on viscous 3-1-1-dimensional hydrodynamics with the Glauber initial condition. Based on our 
results, we argue that the correlation from resonance decays, formed at a late stage of the evolution, 
produce significant effects. In particular, their contribution to the event averages of the coefficients 
of the expansion in the Legendre basis explain 60-70% of the experimental values. We have proposed 
an accurate way to compute these coefficients, independent of the binning in pseudorapidity, and 
tested a double expansion of the two-particle correlation function in the azimuth and pseudorapidity, 
which allows us to investigate the pseudorapidity correlations between harmonics of the collective 
flow. In our model, these quantities are also dominated by non-flow effects from the resonance 
decays. Finally, our method can be used to compute higher-order cumulants for the expansion in 
orthonormal polynomials [1] which offers a suitable way of eliminating the non-flow effects from the 
correlation analyses. 

PACS numbers: 25.75.-q, 25.75Gz, 25.75.Ld 


I. INTRODUCTION 

The mechanism of energy deposition in relativis¬ 
tic nuclear collisions is a subject of intense studies. 
Whereas most of the investigations are concerned with 
the entropy-deposition profile in the transverse plane 
and the resulting transverse expansion (for reviews see, 
e.g., [aEI), the dynamics in the longitudinal direction is 
less explored, and has been recently gaining more atten¬ 
tion with the new experimental analyses from the CERN 
Large Hadron Collider (LHC) expected shortly. Such 
studies could give valuable insight into the initial energy 
and momentum distributions in rapidity [1] , the longitu¬ 
dinal collective dynamics [3, or hydrodynamic fluctua¬ 
tions [S]. Correlations in (pseudo)rapidity can be studied 
in various ways, in particular, as correlations of the trans¬ 
verse flow at different rapidity bins, or as multiplicity 
correlations in rapidity. The first case requires an inter¬ 
mediate collective expansion stage producing flow [ZHS], 
whereas the particle distribution and multiplicity corre¬ 
lations in rapidity are not modified significantly during 
the fireball expansion, thus are expected to reflect more 
closely the initial conditions in the fireball. In other 
words, the multiplicity correlations arise even without 
any collective expansion. 

Correlations of the multiplicity of particles observed 
in high energy collisions in different pseudorapidity in¬ 
tervals have been studied in a number of colliding sys- 
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terns mini El. The most common approach is based 
on the correlation of the number of particles in forward 
and a backward pseudorapidity bins, (npriB), or related 
observables. 


Bzdak and Teaney have proposed to expand the two- 
point correlation function in pseudorapidity in a basis of 
orthogonal polynomials m- The correlations are then 
written in terms of the corresponding expansion coeffi¬ 
cients {anara)- The extracted coefficients can serve to 
parametrize event by event fluctuations of the particle 
distribution in pseudorapidity. A basis of the Legen¬ 
dre polynomials [12] has been used for the expansion of 
the correlation in pseudorapidity for the case of Pb-|-Pb 
collisions at -^/s = 2760 GeV, recently measured by the 
ATLAS Collaboration [T3|. 


In this work we present predictions of the relativistic 
hydrodynamic model for the two-particle correlations in 
pseudorapidity, focusing on correlations generated in the 
late stage of the collision via resonance decays. Our ap¬ 
proach consists of a Glauber Monte Carlo model with 
asymmetric longitudinal emission profile for the initial 
state, and the viscous 3-l-lD hydrodynamic evolution of 
the fireball, followed by statistical hadron emission at 
freeze-out. Our main result is that the late-stage correla¬ 
tions from resonance decays contribute largely (about a 
half of the measured values) to the correlations extracted 
in terms of the (a„am) coefficients. The missing strength 
should be attributed to the correlations generated in the 
earlier stages of the evolution (initial state, jets). 
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II. TWO-PARTICLE CORRELATION 


The two-particle correlation in pseudorapidity, scaled 
by the one-particle distributions, is defined as 


C{vi,V2) 
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where N{ri) denotes the distribution of the number of 
hadrons at 77 and the averaging is over events in a se¬ 
lected centrality class. In the experiment, the correla¬ 
tion function is constructed as the ratio of the histogram 
for particle pairs from physical events to the histogram 
constructed from mixed events in the same centrality 
class [13] 
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We note that the definition ([I I corresponds to the 
scaled second factorial moment of the multiplicity distri¬ 
bution, which depends on the centrality definition and 
the width of the centrality bin. To reduce the effects of 
the overall multiplicity fluctuations, the ATLAS collabo¬ 
ration uses a modified correlation function 
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with Cp{rii) = ^ /Ty C(??i, ? 72 )d'' 72 - Tbe experimental 
analysis suggests that CAr(r7i,772) is approximately inde¬ 
pendent of the definition of centrality [T2j [T3| . 


III. EXPANSION IN ORTHONORMAL 
POLYNOMIALS 


As the shape of the distribution function N{rj)/{N{r])) 
fluctuates event by event, it can be expanded in a basis 
of orthogonal functions [lO] 
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For the case of Legendre polynomials Pn{x), the nor¬ 
malized functions are T„(^) = y^ ^"+^ Pn(a:) [T^ . where 

[—y, Y] is the pseudorapidity range on which the correla¬ 
tion functions are measured, such that the orthonormal¬ 
ity condition takes the form 

J_y (f) (y) y = 


The event-average (andm) can be calculated from the 
two-particle correlation function 



( 6 ) 


The procedure is rather complicated, as first the two- 
particle correlation function must be constructed with 
sufficiently fine binning. In the case of low statistics, 
large binning of C introduces biases. 

The estimate of the integral Q can be simply obtained 
from 
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where the sums are over hadrons in the given event and 
the averages are over events. Equation Q produces very 
stable results, free of the binning bias. 

In the experimental analysis of Ref. [T3] the function 
Cat( 771 , 772 ) instead of C{T]i,r] 2 ) is used in Eq. Q. We 
have checked that in our case the resulting difference for 
the (a„am) coefficients for 1 < tt,, ttt < 5 is very small, 
a fraction of percent Q hence in the following we will use 
( 7 ( 771 , 772 ) in Eq. 0. In addition, the function ( 7 ( 771 , 772 ) 
is, in the experiment, normalized to 1. To conform to 
this convention we rescale the coefficients obtained from 
Eq. 0: 
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In practice, for centrality bins in the model calculation 
defined by the number of participant nucleons, the cor¬ 
rection to the normalization is less than 2%. 

The motivation of the studies of Ref. nnmana was 
to transform the two-particle distributions into a series 
of coefficients (a„am) with a simple interpretation. For 
instance, the coefficient (oioi) is related to the asymme¬ 
try in the entropy deposition in rapidity from the forward 
and backward going participant nucleons. The asymme¬ 
try of the deposition in rapidity is visible in the charged 
particle distribution in pseudorapidity in asymmetric col¬ 
lisions m and in forward-backward multiplicity distri¬ 
butions m- 


IV. RESULTS FROM THE HYDRODYNAMIC 
MODEL 

We use the 3+1-dimensional viscous hydrodynam¬ 
ics m to model the evolution of the fireball created in 
Pb-|-Pb collisions at y/s = 2.76 TeV. The initial entropy 


A correction, which is tiny, could be worked out along the lines 
of Ref. | 12| . 
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C{t 1^, Tig) c=60-70% 



FIG. 1. (color online) Two dimensional correlation function 
in pseudorapidity for charged particles in Pb+Pb collisions at 
y/s = 2.76 TeV at three centrality classes 0-5%, 30-40% and 
60-70% in panels (a), (b) and (c), respectively. 


density in the transverse plane is calculated in GLIS- 
SANDO [ 12 ], implementing the Glauber Monte Carlo 
model. The initial profile in the longitudinal direction 
(in space-time rapidity) is considered in two qualitatively 
different scenarios. In the first scenario, the entropy dis¬ 
tribution in space-time rapidity from a left- and right¬ 


going participant nucleons is of the form 
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and ?/beam is the rapidity of the beam. For the LHC, 
the parameters determining the shape are cr = 1.4 and 
rjp = 2 . 4 . The asymmetric distribution of the de¬ 
posited entropy between the forward and backward ra¬ 
pidity hemisphere leads, together with the fluctuations 
in the number of participants, lead to nontrivial correla¬ 
tions between forward and backward rapidity bins, both 
in multiplicity m and in the flow angle orientation mm- 
The latter has been termed the torque effect, hence we 
label our calculations based on Eq. ( 11 ) as torque. 


The reference scenario assumes that the initial entropy 
profile in space-time rapidity is symmetric. 


/±(? 7 ||) = ^(??||) for |?7|||<2/b earn ■ (11) 

In that case (labeled no torque) in each event the fire¬ 
ball density has a backward-forward symmetry, hence no 
shape fluctuations of odd reflection symmetry are pos¬ 
sible. To summarize, the torque case includes certain 
initial-state fluctuations in rapidity, while the no-torque 
case does not. 

At freeze-out, hadrons are emitted, and later reso¬ 
nance decays occur. The decays of resonances introduce 
short-range correlations of length of about one unit in 
pseudorapidity, leading to a nontrivial structure of the 
two-dimensional correlation functions. Note that another 
source of correlation in the late stage, unrelated to the 
fireball shape fluctuations, is due to local charge conser¬ 
vation [Illlll]. 

The correlation function is calculated from real¬ 
istic, finite-multiplicity events, generated after the hy¬ 
drodynamic evolution with THERMINATOR [T 9 |. We 
use the freeze-out temperature T/ = 150 MeV. The sim¬ 
ulated events include the short-range correlations from 
resonance decays. In Eig.j^we show the two-dimensional 
correlations for three different centrality classes. Gharged 
particles with pj_ > 0.5 GeV and \r]\ < 2.5 are taken to 
simulate the ATLAS acceptance. 

For C{r]i,r]2), plotted in Fig.[^for three sample central¬ 
ities, a peak from the short range correlations is clearly 
visible around ryi ~ 772. When passing to ^^7(771,772), we 
note that the denominator in Eq. ([^ is smaller than one 
at large 1 771,2!, hence it causes relative enhancement of the 
correlation measure in this region. As a result, the shape 
of the correlation function is changed significantly when 
passing from (7(771,772) to Cn{ rji, 112), cf. Eigs.j^an di In 
particular, a saddle-like form appears, corresponding to 
a term of the form ^(77^ — 77^), where ri± = 771 ±772. Such 
a term is expected from event-by-event asymmetry of the 
initial distribution function nni, giving a nonzero value of 
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CN(T1^.'n2) C=0-5% 



Cn{\ti^) c=30-40% 



Cn{\ti^) c=60-70% 



FIG. 2. (color online) Same as Fig. but for the corrected 
correlations function C'jv(? 7 i,? 72 ) ©• 
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FIG. 3. (color online) Calculated coefficients \J{a„an) for 
centrality 0 — 5% (panel a) and 30 — 40% (panel b) for the 
torque and no-torque models, as well as for the oversampled 
events for the torque case (see text for details). 
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predicted by the model reaches about 60 - 70 % of the val¬ 
ues observed experimentally m- The trend of the de¬ 
pendence on the rank n is similar as in the experiment. 
Similar conclusions can be made for the non-diagonal co¬ 
efficients i/—(o„a„_|_2) shown in Fig. W With the avail¬ 
able statistics, we cannot calculate hi^er order averages 
{andn+i) ■ The results for the two scenarios of the the ini¬ 
tial conditions, torque and no-torque, are shown. Inter¬ 
estingly, both calculations give very similar results. This 
shows that in our model the dominant contribution in the 
observed signal comes from the short-range correlations 
due to resonance decays. 


(oiOi). Without this asymmetry, the short-range corre¬ 
lations are expected to be a function of 1 771 —772 1 only [ 5 D]. 

However, in our simulations almost the same value 
of (oioi) is obtained from the correlation functions 
C{r]i,ri2) and C'jv( 7 ?i, 7 ? 2 )- It thus suggests that the ob¬ 
served dependence of the correlation function on ri± does 
not directly prove the existence of correlations induced 
by the event-by-event fluctuations of the distribution. 

In Fig. N we show the calculated coefficients \/JcLtiO^, 
n = 1 ,.. .77 for two sample centralities. The magnitude 


V. DOUBLE EXPANSION OF CORRELATIONS 
FUNCTIONS IN AZIMUTHAL ANGLE AND 
PSEUDORAPIDITY 

Collective flow in ultrarelativistic heavy-ion collisions 
causes all particles to be emitted in a correlated way, 
which leads to azimuthal asymmetry in hadron distribu¬ 
tions. The correlation function between two pseudora¬ 
pidity bins, constructed for multiplicity correlations as 






































5 


0.010 

0.008 

: 0.006 

« 

0.004 

I 

^ 0.002 
0.000 


0.014 
0.012 
- 0.010 
J 0.008 
«= 0.006 
^ 0.004 
^ 0.002 
0.000 

1 2 3 4 5 


-. 

. 

. 

(a) 

■ TORQUE 
□ NO TORQUE 

• OVERSAMPLED 

■ 


c = 0-5% 







• 

• 

• 

• 

• 

1 

2 3 

4 5 

(b) 

hydro 3+ID 

Pb +Pb @2.76 TeV 

% 

■ . 

c = 30-40% 


0 

% 

■ 


0 



■ 

• 

• • 

• • 





I 



I 






U.UU8 

(a) 


■ 

TORQUE 


0.006 


□ 

NO TORQUE 




• 

OVERSAMPLED 


□ 



c=0-5% 


0.004 


n 

m 



0.002 

• 



m 

□ 

0.000 


• 

• 

• • • 


0.012 

1 

2 

3 

4 5 6 

7 

-1 ■ ■ • 




1 

(b) 


hydro 3+lD 


0.010 


Pb 

+ Pb@ 2.76 TeV 


0.008 

■ 

: □ 

□ 

■ 


c =30-40% 

0.006 

; 


% 



0.004 

• 



e % 






£ 

i 

0.002 






0.000 


• 

• 

• • • 

• 


1 

2 

3 

4 5 6 

7 


n 


n 


FIG. 4. (color online) Same as in Fig. bnt for the the 
coefficients (a„a„+ 2 ). 


FIG. 5. (color online) Same as in Fig. but for the the 
coefficients of the second order harmonic ^ [an[2\ a„[—2]). 


in the preceding sections, can be straightforwardly gen¬ 
eralized for each harmonic flow component for any two 
pseudorapidity bins. Various techniques are applicable 
here. The rapidity dependence could be decomposed into 
principal components | 21 ] . but in practice the principal 
component analysis may be difficult and restricted to the 
lowest eigenmodes. Alternatively, the harmonic flow cor¬ 
relations in pseudorapidity can be expanded in a basis 
of suitable orthogonal polynomials, in full analogy to the 
multiplicity case. This provides insight into a different 
characteristic of the flow, with possibly different sensi¬ 
tivity to non-flow effects than in the multiplicity correla¬ 
tions discussed in Sec. El 

Let us define the correlation coefficients for the n-th 
order harmonic flow as 


{aj[n]ak [-n]) 


T, (^) Tk if) \ 


( 12 ) 

In the above equation, use we the normalization of the 
correlations function by l/(A^(? 7 i))(V( 772 )) as in Eq. ([^, 
but the formula can be written analogously for the corre¬ 
lation function of flow vectors in two rapidity intervals as 
used in m- Note that the linear part of the pseudora¬ 


pidity dependence of the torque effect for the orientation 
of the flow angle [ 7 ] contributes to the (ai[n]ai[—n]) co¬ 
efficient. 

In Figs. and we show the decomposition coeffi¬ 
cients © of the elliptic and triangular flow correlation 
at different pseudorapidities. We compare calculations 
using the torque and no-torque scenarios for the initial 
conditions, as in Sect. |IV| We notice that the two calcu¬ 
lation give similar results, although in the no-torque case 
the odd coefhcients should vanish within the statistical 
uncertainties. Our results mean that in the decomposi¬ 
tion of the flow correlations in pseudorapidity, the dom¬ 
inant contribution comes from resonance decays. The 
same effect has been noticed in the analysis of factor¬ 
ization breaking for flow at different pseudorapidities [ 5 ] 
(the torque effect). 


VI. HIGHER-ORDER CUMULANTS 

Non-flow correlations have a significant contribution 
to the measured (a„am) coefhcients. In this Section we 
show results of an idealized calculation with the non-flow 
effects removed. The coefficients are calculated using 
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FIG. 6. (color online) Sa me as in Fig, [^ but for the coefficients 
of the third harmonic -y/(a„[3]a„[—3]). 


oversampled events, where for each hydrodynamic evo¬ 
lution several hundreds of THERMINATOR events are 
generated and combined together. That way non-flow 
effects are damped. The procedure is equivalent to a 
Monte Carlo integration of one-particle densities in each 
event. The results for (aiOi) and ( 0202 ) are show in 
Fig.0 The genuine effect due to event by event fluctu¬ 
ations of the rapidity distributions is small. The coeffi¬ 
cients from the shape fluctuations are much smaller than 
correlations from resonance decays. The dependence on 
the rank n of the coefficients {anttn) from oversampled 
events is presented in Figs. [^to[^ As expected in the 
torque scenario, involving forward-backward asymmetry, 
(aiOi) and (ai[n]ai[—n]) have the largest magnitude. 


Two particle correlations from resonances are removed 
in higher order cumulants [1], while genuine correlations 
due to initial-state fluctuations of rapidity distributions 
0 do contribute. There are many possible combina¬ 
tions of the fourth-order cumulants that can be used for 
the purpose. In particular, one can define the simplest 
fourth-order cumulant for the multiplicity fluctuations 
as [5 


(a|)c 


e' 
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(13) 


where the subscript c stands for the connected part and the prime denotes summation over different particles. For 
the flow correlations in pseudorapidity, the most general cumulant is of the form 


/ _ T ( '!hx\ pi'mi'pa, T- ( I 
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with X]fc=i = 0- The simplest fourth-order cumulants are 
{ak[n]ak[n]ak[-n]ak[-n])c = 
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We have attempted to compute the fourth-order cumu¬ 
lants (13) and (|15|) in our simulation, however, with the 


available statistics (20000 events in each centrality class) 
the statistical errors are too large, of the same order as 
the square of the second order cumulant. The application 


of the cumulant method [T] is possible on large-statistics 
experimental data and the results could be compared to 
model calculations using one-particle densities (such as 
the results for the oversampled events presented above) 
that neglect the non-flow effects. 
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FIG. 7. (color online) Coefficients \J (aioi) (panel a) and 
( 0202 ) (panel b) for the torque model and for the oversam¬ 
pled events for the torque case, plotted as functions of the 
number of participants (determining centrality). The ATLAS 
data come from Ref. |13 |. 


that the correlation from the resonance decays, formed 
at a late stage of the evolution, produce significant ef¬ 
fects. In particular, their contribution to the coefficients 
(onOm) in the expansion of the correlation function in the 
Legendre basis take 60-70% of the experimental values. 

While our model incorporates only some possible 
sources of correlations (those from the torque effect and 
the resonance decays), it shows their relevance in the 
analyses. Other, not incorporated non-flow effects, in¬ 
clude the hadron production from jets, local current con¬ 
servation, or additional sources of rapidity fluctuations 
in the initial state [5] . 

On the methodological level, we have proposed a new 
way to compute the (a„am) coefficients, independent of 
the binning in pseudorapidity, and applied it in our sim¬ 
ulations. Also, we have developed a double expansion of 
the correlation function in the azimuth and pseudorapid¬ 
ity, which allows to probe and quantify the rapidity cor¬ 
relations between harmonics of the collective flow. We 
have found that in our model these quantities are also 
dominated by non-flow effects. Our method can be used 
for higher-order averages of the orthogonal polynomials, 
in particular for cumulants. This offers a way of eliminat¬ 
ing the non-flow effects [1] , but requires very large statis¬ 
tics, which, fortunately, is available in the experiments. 
These measures could be compared to model calculations 
with oversampled events, where sufficient statistics can 
be achieved. 

We note that a study using similar methods and lead¬ 
ing to similar results has been independently and simul¬ 
taneously presented in Ref. [22] . 
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